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We demonstrate a mesoscopic self-collimation effect in photonic crystal superlattices consisting 
of a periodic set of all-positive index 2D photonic crystal and homogeneous layers. We develop an 
electromagnetic theory showing that diffraction-free beams are observed when the curvature of the 
optical dispersion relation is properly compensated for. This approach allows to combine slow light 
regime together with self-collimation in photonic crystal superlattices presenting an extremely low 
filling ratio in air. 
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Electrodynamic theory in negative-index materials, 
originally introduced by V.G. Veselago in 1960 pLi], has 
known impressive progress since the advent of metama- 
terials and photonic crystals [2]. Both approaches en- 
able an invaluable light control which originates for the 
former from the realization of effective electric and mag- 
netic permittivities and for the latter from versatile pho- 
tonic dispersion relations. Many concepts developed for 
metamaterials in linear and nonlinear regimes have been 
transposed to PhCs devices since they enable the design 
of devices working in the visible range of frequency [3|-|5| . 
The frontier between PhCs and metamaterials has be- 
come fuzzier with recent emergence of Bragg mirrors in- 
cluding negative-index materials. These composite one- 
dimensional PhCs, alternating positive- and negative- 
index layers, have in particular shown an intriguing pho- 
tonic band gap when around a particular frequency, the 
average refractive index over one lattice period is null [6- 
fTlj . This forbidden range of frequency named zero-n gap 
is for example insensitive to the thickness of the lattice 
period or disorder pjj] and electromagnetic signals can 
be tunneled without any phase delay because of a phase 
compensation mechanism [Ts*]. Resonant modes lying in 
the zero-n gap can furthermore appear without breaking 
the lattice symmetry when a Fabry-Perot optical con- 
dition is satisfied [6^. Surprisingly, it has been demon- 
strated that zero-average index metamaterials support- 
ing these resonant modes can either propel self-collimated 
beams or focalize optical signals [14]. These properties 
have been explained by introducing a harmonic average 
index parameter that takes null or negative values while 
the average index is kept to zero. However, this approach 
may be difficult to implement experimentally since an ac- 
curate control of the metamaterial optical dispersion is 
required [14, 15]. This drawback is circumvent in an al- 
ternative device proposed by V. Mocella and co-authors 
which consists of a periodic stack of 2D photonic crystals 
(PhCs) and air-layers of equal thicknesses pJi] . The PhC 



slabs were designed to behave as fiat lenses of -1 effective 
index originating from an isotropic dispersion relation 
obtained with the second photonic band. In that case, 
self-collimation can be interpreted as successive focaliza- 
tions of the beam by a lensing effect. Experimental data 
also reveal that sub-wavelength beams can efficiently be 
relayed over millimeter-scale distance thanks to a smart 
optimization of the PhC interfaces. This quasi-perfect 
impedance matching prevents however to open a large 
zero-n gap so that self-colHmation was observed close to 
the band gap edge. 

Nevertheless, to date, this beam shaping mechanisms 
is restricted to photonic devices realizing a quasi-zero- 
average index condition by the use of negative-index ma- 
terials and this is achieved in structures with high filling 
factor in air (of 76%) [16^. In this Letter, we demon- 
strate that mesoscopic self-collimation also arises in pho- 
tonic crystal superlattices of all-positive index materials 
and when the zero-average index condition does not hold. 
Diffraction-free beams are demonstrated when the curva- 
tures of the photonic dispersion curves of the successive 
media are offset. We propose an electromagnetic the- 
ory that enables the design of Bragg media of extremely 
low filling factor in air or presenting slow self-collimated 
light. Let us start with an infinite 2D PhC consisting 
of a square lattice of air holes (lattice constant a, radius 
r/a = 0.2) etched in a dielectric medium of optical index 
2.9. This 2D PhC has an air filling factor of around 12%. 
The first photonic band, computed with the plane wave 
expansion method [18], shows fiat iso- frequency curves 
(IPCs) of quasi null curvature around the reduce fre- 
quency a/ A = 0.235, Fig.l. At this frequency, all the 
waves constituting a beam propel with parallel group ve- 
locities (vg = gradkCo'(k)) pointing in the TM direction 
revealing the self-collimation effect [17]. At lower fre- 
quencies (red area in Fig.l), the IPCs have a positive 
curvature leading to a diverging beam, as in an homo- 
geneous medium. Conversely, for higher frequencies, the 
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FIG. 1. (color). Iso- frequency curves of the first band for TM 
polarization. The arrows indicate the direction of Vg and the 
red and blue areas correspond respectively to IFCs of positive 
and negative curvatures. 



negative curvature results in beam focusing. Note that in 
this regime, negative refraction is not linked to a negative 
phase index related to opposite phase and group veloci- 
ties since here k.Vg > for all Bloch waves la. These 
remarks show the crucial role of the local curvature of 
the IFCs for beam shaping operations [20]. 

Suppose now a PhC superlattice consisting of a peri- 
odic set of the previous 2D PhC sized in slabs of thickness 
di and of homogeneous layers of thickness d2 and refrac- 
tion index n2 = 2.9. The TM direction of the 2D PhC is 
chosen parallel to the direction of the PhC superlattice 
stack. Mesoscopic self-collimation should appear when 
the positive and negative curvatures of the dispersion re- 
lations of the homogeneous and the PhC layers are ex- 
actly balanced. At the scale of the Bragg crystal macro- 
period, D = di -\- the resulting average curvature 
should vanishes provided that appropriate layers thick- 
nesses are chosen. This concept is developed in the frame- 
work of a beam propagation model. Consider a Gaussian 
beam of initial waist Wo launched in a PhC superlat- 
tice presenting a ratio in PhC layers a = d2/di. After a 
unit cell distance the beam is derived by the inverse 
Fourier transform of the product of the initial Gaussian 
beam envelop U'{k^) = Wo/(20r) exp(-(A:^Wo/2)^) by 
the phase propagators of each media: 

U{x,D) = TF-^ {U'{k^)Pphc{k^,di)Phom{kx,d2)} 

(1) 

where the general expression of the phase propagator 
in a medium of dispersion relation ky{kx) is P{kx^y) = 
e-x.Y)[iky{kx)y\- In high symmetry directions (such as in 
VM direction), the parity of the IFCs imposes a null first 
derivative dky/dkx- Consequently, the Taylor expansion 
of the wave vector k,, written as: 



ky{kx) — ky{kxo) {kx kxo) 



2k 



(2) 



The right hand term is linked to the local curvature of 



the reciprocal basis is chosen as kxo = in TM direction 
and we define kyo = ky{kxo)' In that case, for the ho- 
mogeneous slab, the curvature is related to the refraction 
index by 1/khom = 1/(^0^2) (^0 being the wavenumber 
in vacuum). Similarly, let introduce for 2D PhC struc- 
tures the dimensionless parameter ric named curvature 
index which is proportional to the curvature radius of 
the IFCs: Uc = kphc/ko. As depicted in Fig. 2a, l/uc 
vanishes at the self-collimation frequency a/X = 0.235 
and takes positive or negative values at lower or higher 
frequencies respectively. The expansion of Eq.(j2j) (assim- 
ilated to the paraxial assumption) allows one to express 
the phase propagator over one unit cell: 



Pceiiikx.D) = exp 



(3) 



which is simply the product PhomPphC- This phase prop- 
agator involves the average phase {kyo) and the harmonic 

average curvature index given by ^7^^ = ^nJT)' 
Using Eq.([T]) and Eq.(|3]) one can find an analytical ex- 
pression of the beam inside the Bragg media. The waist 
of the Gaussian beam after one macro-period D depends 
in particular on the average curvature: 



W{D) = Wq\1- 



OqDV / 1 ' ^ 



Wo 



(4) 



where Oq = A/(7rWo) is the divergence angle [2T|. This 
result demonstrates that the initial beam waist is re- 
covered after propagating through one Bragg period, 
W{D) = Wo, when the average curvature is zero: 



di do 

— + — = 

Tie n2 



(5) 



the IFC by i = 



d^ky 
dkx^ 



For the sake of simplicity 



Note that here mesoscopic self-collimation is based on 
the exact balance of the local IFC curvatures of each 
layers and does not req uire to open a zero-n gap that 
implies (kyo) = [1, [ij. The beam propagation model 
in PhC superlattices provides a simple method to deter- 
mine the self-collimation frequency from the 2D infinite 
PhC dispersion properties. The working frequency is in- 
deed fixed by the intersection point of the local curvature 
l/uc with the straight line —a/n2 deduced from Eq.([5]), 
Fig. 2. For ratios in PhCs layers a = 1 and a = 3, meso- 
scopic self-collimation should respectively be obtained at 
the reduced frequencies a/X = 0.241 and a/X = 0.244, 
Fig. 2a. The operating frequency is indeed blue shifted 
when increasing a because the PhC layers must focus 
more strongly to cancel diffraction in thick homogeneous 
slabs. It is also worth noting that these frequencies are 
invariant with respect to the macro-period D provided 
that Eq.(|5|) is satisfied. This invariant scaling property 
has been checked by calculating the photonic band dia- 
grams of PhC superlattices with the plane wave expan- 
sion method. Flat IFCs are found for a constant a and for 
several di varying from aV2 to 4a>/2 (by step of a>/2). 
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FIG. 2. Curvature function l/ric computed in VM direc- 
tion versus the reduce frequency, (a) For the infinite 2D 
PhC described in Fig.l and of unit cell depicted in the in- 
set. The intersection points between l/ric and the straight 
dashed- lines for a = 0, 1 and 3 correspond respectively to the 
self-collimation frequencies for the unbounded 2D PhC and 
for PhC super lattices of two filling ratios in PhC layers, (b) 
Graph of l/ric computed for the supercell corresponding to a 
PhC superlattice with a — 1 and di — 3a\/2 shown in the in- 
set. Mesoscopic self-collimation arises for a/A = 0.239 where 
l/uc vanishes. 



For example, it is seen in Fig. 2b that for di = 3aV^ 
and a = 1, the curve l/ric vanishes at a/ A = 0.239 
in agreement with the beam propagation model. The 
weak difference (smaller than 1%) between this operat- 
ing frequency and the prediction of the propagation beam 
model shows that our model is thus accurate enough to 
guide designing and fabrication of self-collimated PhC 
structures. To corroborate these results, the propaga- 
tion of an incident beam has been performed with 2D 
FDTD simulations through PhC superlattices presenting 
two filling ratios a [23| . We found a self-collimated beam 
over a length of 400a at the reduce frequencies 0.239 and 
0.243 for a = 1 and 3 respectively. Fig. 3. Theses results 
demonstrate that mesoscopic self-collimation is achieved 
in layered structures of extremely low filling factor in air. 
The air filling factor given by / = 7r(r/a)^/(l + a) in the 
superlattices of Fig. 3a and 3b is respectively divided by 2 
and 4 versus to the initial unbounded 2D PhC. Compared 
to the bulk semiconductor case where the beam waist at- 
tains 3 times Wo after the same distance 400a (Fig. 3c), 
the beam keeps the same initial waist for the PhC super- 



FIG. 3. (color) Map of the magnetic field modulus for a Gaus- 
sian beam of waist Wo = 8a (a) At the reduce frequency 
a/\ — 0.239 for a PhC superlattice for di — 4a^/2 and a = 1. 

(b) At the reduce frequency 0.243 for di — 3a>/2 and a — ?>. 

(c) The same Gaussian beam propagating in a homogeneous 
medium of optical index n2 = 2.9 and for a/A = 0.239. 



lattice presenting an air filling ratio of only 3% (Fig. 3b). 
This value could furthermore be decreased by considering 
2D PhC layers of lower curvature index. The selectivity 
of PhC superlattices can be evaluated by computing the 
range of frequency supporting the self-collimation effect 
in the same way than in [20]. Here, in the framework of 
the beam propagation model, it is related to the deriva- 
tive of the waist versus the reduce frequency which is 
derived for a large propagation distance D from Eq. (j4]) 
: dW/d{a/\) = Oodid{l/nc)/d{a/ X). The selectivity is 
then driven by the slope of the curvature l/ric of the 
infinite PhC and it increases with higher values of a/ A, 
see Fig. 2a. For example, the self-collimation bandwidth 
is two times smaller for a = 1 compared to the infinite 
PhC case. Conversely to the approach proposed in |2Q| , 
the self-collimation bandwidth is reduced for PhC super- 
lattices working with the first photonic band. This con- 
clusion may however be reconsidered for others photonic 
dispersion relations. 

Besides of these linear properties, these PhC superlat- 
tices present also attractive properties for nonlinear ap- 
plications such as laser emission or frequency conversion. 
Nonlinear optical processes developed for Bragg mirrors 
can be revisited owing to the original optical properties of 
PhC superlattices. First, the latter enables to guide self- 
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FIG. 4. (color), (a) Band diagram in the FM direction of the 
PhC superlattice for a = 1 and di — depicted in the in- 

set. The red circles indicate the location of the self-collimation 
frequency 0.239. (b) Group index versus the reduce frequency. 



collimated beam in devices presenting a strong amount of 
non-etched semiconductor material. Indeed, as compared 
to structures reported to date [l6| , the air filling factor is 
reduced by a factor 25 resulting in a potentially 97% of 
optically active or nonlinear semiconductor. Moreover, 
as seen in Fig. 3b, the PhC superlattice of low ratio in 
PhC layers behaves like a set of coupled cavities where 
a strong electromagnetic field localizes. Second, light- 
matter interaction is also enhanced when light is slowed 
down at frequencies close to a photonic band gap (PBG) 
edge [22]. In lasers, the optical gain increases linearly 
with the group index, Ug = c/vg, and the efficiencies of 
second or third harmonic generation nonlinear processes 
are enhanced as n^. Unfortunately, self-collimation in 
2D or 3D PhCs arises far away from the PBG edge pre- 
venting thus to work in a slow light regime. In our case, 
the infinite square lattice of air holes presents a complete 
PBG between the first and the second bands spanning the 
frequency range 0.247-0251. Self-collimation arising at 
the reduce frequency 0.235 is associated to Bloch modes 



of around 2.96 group index which is almost the refrac- 
tion index of the host material. To increase this low 
group index one can play with the additional degree of 
freedoms provided by PhC superlattices. Since these lay- 
ered media are basically Bragg mirror structures, we can 
open an additional PBG at the vicinity of the mesoscopic 
self-collimation frequency by an appropriate choice of the 
macro-period D and the parameter a. For example, the 
band diagram computed in the TM direction, for a = 1 
and di = 4a\/2 corresponding the structure of Fig. 3a, re- 
veals that a new PBG centered at a/A = 0.237 and orig- 
inating from Bragg scattering opens. Fig. 4a. With these 
parameters, the self-collimation frequency, a/A = 0.239, 
appears now at the band gap edge located at the reduce 
frequency 0.2385. A high group index of around 50 is ob- 
served for a PhC superlattice of 6% air filling factor that 
combines both slow light and the self-collimation effects, 
Fig.4b. 

In summary, a theory of photonic dispersion cur- 
vature compensation has been presented to predict 
self-collimation of light in PhC superlattices. The 
beam propagation model is shown to be relevant for 
the conception of PhC superlattices and gives an 
insight on the optical mechanism behind the meso- 
scopic self-collimation effect. This theory reveals that 
negative-index materials are not necessary to generate 
diffraction-free beams leading thus to a general approach 
for the design of Bragg mirrors based on 2D photonic 
crystals properties. Mesoscopic self-collimation is 
demonstrated in PhC superlattices presenting 97% of 
non-etched semiconductor. Slow self-collimated light 
is in addition showed in appropriately designed PhC 
superlattices that work at a frequency located at the 
edge of an additional photonic band gap. To our opinion, 
these novel properties pave attractive avenues for the 
design of active or nonlinear devices presenting beam 
shaping functionalities. 
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